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Abstract. We relate the pressure 'dip' observed at the bottom of a sandpile prepared by successive 
avalanches to the stress profile obtained on sheared granular layers in response to a localized vertical 
overload. We show that, within a simple anisotropic elastic analysis, the skewness and the tilt of the re- 
sponse profile caused by shearing provide a qualitative agreement with the sandpile dip effect. We conclude 
that the texture anisotropy produced by the avalanches is in essence similar to that induced by a simple 
shearing - albeit tilted by the angle of repose of the pile. This work also shows that this response function 
technique could be very well adapted to probe the texture of static granular packing. 

PACS. 45.70.-n Granular systems - 45.70.Cc Static sandpiles - 46.25.-y Static elasticity 



The stress distribution below a pile of sand has been 
one of the problematic issues of the statics of granular 
materials in physics over the last few years P^. In fact, 
experiments have shown that, when a granular pile is pre- 
pared from a point source, the bottom pressure profile has 
a clear local minimum - a 'dip' - below the apex [21314) . 
The existence of this pressure dip has been strongly de- 
bated, and it is now well established that the presence or 
absence of this dip is closely related to the preparation 
history of the pile. This was demonstrated by Vanel et 
al. U. Using the same sand and experimental apparatus, 
these authors could generate the stress dip using a local- 
ized deposition technique or cause the dip to disappear 
by constructing a pile in successive horizontal layers. Sim- 
ilar conclusions were reached for two dimensional heaps 
with photo-elastic grains |S] , and in numerical simulations 

This interesting effect has inspired the development 
of new models to describe how forces are transmitted in 
dense granular materials. Among them are those proposed 
by Bouchaud et al., initially developed in the context of 
the sand pile dip [911011 1112"] . and further extended to 
other geometries like that of the silo |12I13| . This approach 
is also intended to describe a collection of systems includ- 
ing dense colloids, granular matter or foams ^J. At the 
macroscopic level, these features are modelled by hyper- 
bolic, partial differential equations (PDE) for the stress 
tensor. Although no explicit link was established, the char- 
acteristics of these hyperbolic equations were intuitively 
thought to be related to the mesoscopic 'force chain' net- 
work whose structure and orientation were shaped by the 
previous history of the granular assembly - see also |14I15| 



concerning force chains. Plasticity theories for granular de- 
formations are also of hyperbolic type, although concep- 
tually different than the previous cited models. From the 
classical, soil mechanics point of view, below the plastic 
threshold, granular material is thought to behave as an 
effective elastic material with PDE's that belong to the 
elliptic class |TJ|. Finally, sound wave propagation tech- 
niques and numerical simulations of confined granular as- 
semblies indicate that that assesment of effective elastic 
constitutive relations is still an open and difficult issue 

na- 

In order to distinguish between the very different math- 
ematics of hyperbolic and elliptic PDEs a stress response 
experiment was proposed |9I19[ . For instance, the pressure 
profile measured at the bottom of an isotropic elastic hor- 
izontal layer in response to a localized vertical overload at 
its top surface should be a single broad peak, while hyper- 
bolic models would yield in 2D two thin peaks or a ring in 
3D. Seve ral experiments [^112^12^12^,^2?)] and sim- 
ulations |27I28I29I3UI31I32I33I34|5515d1 have recently ad- 
dressed the issue of the stress response in a granular layer. 
Collectively, they have demonstrated two key points: 
(i) the shape of the pressure response profiles is generally 
not in agreement with the predictions of the hyperbolic 
models. For the generic case of a disordered packing of 
frictional grains, the measured profiles show an elliptic- 
like behavior with one single peak broadening in propor- 
tion to the thickness/depth of the layer - except in |2U[ . 
For well ordered packing, the measured response functions 
may produce rays diverging from the source, but this type 
of response is also compatible with anisotropic elasticity 
13013713^ : 
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Fig. 1. Left: Sketch of the 2D shear box. The system is strained up to an angle <j>. Middle: Visualization of the force chains 
after shearing. The analysis of their orientations shows that a direction at r = 45° is prefered [25) . Right: Average response to 
a vertical overlaod after shearing with a strain angle cj> — 5° . The stress maximum is oriented at (3 — 22° to the vertical. 



(ii) the geometry of the response function is in fact sim- 
pler than that of the pile, and offers rich possibilities. Re- 
sponse function measurements directly detect any symme- 
try breaking due to texture originating from either a spe- 
cific preparation or from any other external action on the 
pile. The layer can be loose, dense, compacted, sheared, 
avalanched, sedimented, ordered, vibrated, and so on. 
Nevertheless, there are still a number of unresolved ques- 
tions. For instance in the limiting case of isostatic pil- 
ings, numerical simulations |39I40) and theoretical argu- 
ments [41142143144115) indicate that a hyperbolic equa- 
tions should describe the stress propagation, although this 
is not in agreement with the work of Roux |45j . 



Experiments have shown that the response to a point 
force is very sensitive to the preparation of the system pH] . 
For instance, for a 2D system which has been subjected to 
strain in a 2D shear box, the response function is skewed 
in the direction of the shear showing a strain-induced 
anisotropy. Here, we extend this result to 3D granular as- 
semblies and propose a relation to the pressure 'dip' ob- 
served at the bottom of a sand pile when prepared by 
successive avalanches. Interestingly, anisotropy induced by 
preparation was suggested by Savage to explain this 
phenomenon. Note that the occurence of a stress solution 
with a dip can also be produced in a model pile composed 
of an elastic core and plastic wings |47j . 



The paper is organized as follows. We first describe ex- 
perimental results obtained for the response of a sheared 
granular layer in 2D and in 3D (sectional and for a layer 
prepared by successive avalanches. We then present a the- 
oretical anisotropic elastic analysis of these experimental 
data in 2D and a numerical analysis in 3D (section 0) . 
Thereafter, we numerically solve the case of a conical heap 
(section |3J). We close with conclusions and suggested per- 
spectives. 



1 Experimental response function on sheared 
granular layers 

The response function experiments of interest here, were 
carried out in shear cell geometries in two and three di- 
mensions. In both cases, the cell consisted of vertical bound- 
aries that could be tilted quasi-statically by an angle <fi 
with respect to the vertical axis. This process deformed 
the samples from an original rectangular geometry to that 
of a parallelogram. 

In 2D, the spacing between the horizontal boundaries 
was maintained strictly constant, so that the sample vol- 
ume also remained constant. The particles were pentagons 
made of a photo-elastic material, which allows a direct 
measurement of the local force response and of the force 
chains. More details on these experiments can be found in 
Geng et al. As we can see from figure^ the salient fea- 
tures of these series of experiments are (i) the force chain 
network is oriented at 45° from the horizontal axis, which 
is the principal compressive direction and (ii) the response 
to a vertical force is tilted with an angle (3 with respect of 
the vertical direction, which illustrates perfectly the sym- 
metry breaking due to shearing. In a previous contribution 
we argue that this r = 45° angle is simply related to the 
principal axes of compression and dilation (respectively 
directions 1 and 2 on figure HJ. 

In 3D, we carry out two experiments to extend the 
study of the influence of shear on the texture of a granular 
assembly. First, we built an apparatus similar to the shear 
cell already used in 2D. Second, we deposited the material 
in a horizontal layer by the superposition of successive 
avalanches. 

The shear cell - The shear cell is sketched in the left 
part of figure In order to probe the mechanical proper- 
ties and the symmetries associated with the induced struc- 
ture, we mounted a capacitive stress probe at the bottom 
of the cell (labeled 4 in figure |2J) ; the applied force on 
the top (labeled 1) is moved horizontally. Specifically, the 
force response measurements use a low frequency modu- 
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Fig. 2. Left: Sketch of the 3D shear box. 1 top lid, 2 tiltable lateral boundaries, 3 fix glass wall, 4 capacitive stress probe. Right: 
Visualization of the grain displacement due to shear. Coloured grains have been put along four vertical lines before shearing. 
These grains are located along inclined lines afterwards, so that the strain field of the layer looks reasonnably homogeneous. 



lation of the localized stress imposed at the top of the 
sample, and lock-in detection of the vertical pressure re- 
sponse on the bottom, as described by Reydellet et al. |21) . 
The lock-in detection provided a large signal/noise ratio, 
which allowed us to apply tiny forces, of the same order 
as the weight of few grains. This prevented any significant 
deformation of the packing which would have modified the 
fragile structure of contacts. A horizontal layer of grains 
(Fontainebleau sand, slightly polydisperse round grains, 
d ~ 300/xm size) was initially prepared by pouring the 
material into shear box of horizontal size I — 25cm. The 
lateral width of the box was w = 20cm. The total height 
h of sand ranged between 4 and 15cm. 

The pouring procedure involved a sieve which was slow- 
ly raised in order to create a uniform rain of sand. There- 
after the sand layer was packed by pushing on the free- 
surface with a plate, and by simultaneously taping on the 
lateral edges of the box. Note that these conditions are 
similar to the so-called 'dense preparation method' pro- 
posed earlier by Serero et al. 22 . A vertical stress re- 
sponse was determined before and after shearing. During 
the shear deformation of the sample, a weight (~ 40kg) 
was imposed on the top surface and the lateral boundaries 
(labeled as 2 on the left of figure [2J were slowly tilted up 
to a final angle (j>. The dense preparation as well as the 
large imposed load seemed necessary to avoid inhomoge- 
neous deformations of the free-surface and help to hinder 
the formation of localized shear-bands. This result was 
confirmed by direct inspection of the displacement fields 
of the grains lying on the transparent side of the lateral 
boundaries. In order to monitor this displacement field, 
three vertical columns of colored grains were inserted next 
to the lateral boundaries prior to shearing (see right panel 
of figure ■ We note that the formation of a shear band 
was hindered only for limited values of aspect ratios (typ- 
ically for ratios h/l from 0.3 to 1 used here) and for shear 
angles smaller than 5°. We also point out that due to the 
rather close-to-unit aspect ratio of the shear box, the re- 
sponse for large horizontal distances between the source 



and the stress probe was not measured so as to avoid the 
influence of lateral walls. 

In figure |3| we show the experimental results for the 
vertical stress response a zz after shear for several shear 
angles 4> and layer depths h. In this figure, the abscisa is 
the horizontal coordinate x normalized by the sand depth 
h. We normalize the stresses a zz (x) by a zz (0) the value 
of the stress at x = (i.e. immediately below the piston). 
This normalization is due to the finite lateral extension 
of the box since, contrary to previous situations |21I22| . 
we cannot obtain a proper normalization of the stress by 
integration (i.e. assess the total force applied). For shear 




rescaled horizontal position x/h 

Fig. 3. Response profiles of a 3D sheared layer of 
Fontainebleau sand for various values of the layer thickness 
h and the strain angle cj>. The maxima of the profiles corre- 
spond to a tilt angle f3 of 8° ± 1°. The collapse of the data is 
resonnably fitted by the SEM (see section |5]for the definition of 
the SEM and the choice of the different parameters, e.g. Poisson 
ratios, ...) with t = 0.67 and u = 2.14 when r is set to 45°. 
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Fig. 4. Sketch of the avalanche preparation procedure. The 
avalanche angle ip = 30°, so that the expected angle of 
anisotropy is r = 15°. 
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-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 
rescaled horizontal position x/h 

Fig. 5. Two response profiles of a 3D layer of Fontainebleau 
sand prepared by successive avalanches. The thickness of the 
layer is h — 8cm. In comparison to figure |3] the shape of the 
profile is also well reproduced by the SEM (idem, see section |5J 
with t — 0.67 and u — 2.14 but now when r is set to 15°. 



angles <p varying from 2° to 5°, the experiments clearly 
show, a skewing of the response function in association 
with a displacement of the maximum corresponding to an 
angle j3 = 8°±1° (deduced from the slope of the empirical 
relation between the height of sand and the shift). Thus, as 
in the 2D case, the response indicates that, due to shear- 
ing, the granular assembly is clearly anisotropic (at least, 
the vertical direction is no longer an axis of symmetry). It 
is worth mentioning that this tendency is robust, as it was 
observed for other types of grains (larger, slightly polydis- 
perse, rough sand grains of diameter ~ 1mm, and smooth 
spherical beads of diameter ~ 1.5mm). 



The avalanche preparation - Next we prepared a layer 
of sand structured by successive avalanches (see figure 0] 
for a sketch of the preparation procedure). The experi- 
mental results for two independent preparations are dis- 
played in figure [S] and rescaled in a way similar to the 
case of the 3D shear box. Again, we observe a skewing of 
the response indicating an anisotropic texture induced by 
the avalanches. However, the tilt effect of the response is 
weaker, as it corresponds to an angle for the locus of the 
maxima j3 — 3.5° ± 1°. 



2 Elastic calculations 

We next turn to modeling the experimental results using 
elasticity theory. We show here that a simple anisotropic 
description is sufficient to capture all the salient features of 
the data. Isotropic elasticity, in both 2D and 3D, has only 
two constitutive parameters for a given material, namely 
the Young modulus E and the Poisson ratio v. When 
anisotropy is introduced, however, the number of inde- 
pendent parameters increases significantly. For example, 
five parameters are required for the case of simple 2D or- 
thotropy that we use below. In 3D, this number is even 
higher. Fortunately, for the static equations, only combi- 
nations of these parameters enter (e.g. the parameters r 
and t, below). Still, a meaningful fit of the response pro- 
files to an anisotropic elastic model is a formidable task. 
Although such an approach is very instructive, it is not an 
essential point for the present purposes. We will present 
the results of this type of approach elsewhere. 



2.1 Elastic response in 2D 

Let us first focus on the simpler two dimensional case. 
What is indicated by the 2D shearing experiment described 
in the previous section, is that the sheared force chain net- 
work clearly has a preferred direction at an angle of 45° 
with respect to the initial axis of the box. Contacts be- 
tween the grains are gained or reinforced in this strong 
direction, whereas they are lost or weakened in the per- 
pendicular direction. It is then expected that an effective 
elastic medium model for this situation would have a stiff 
direction (direction 1) which makes a fixed angle t = 45° 
relative to the vertical axis z (x denotes the horizontal co- 
ordinate) characterized by a Young modulus E\ . Perpen- 
dicularly to direction 1 is the softer direction (direction 2) 
that is characterized by a Young modulus E 2 < E\ . This 
orthotropic elastic description is completed by three addi- 
tional parameters: two Poisson coefficients v\2 and i>2i and 
the shear modulus G. Then, the strain-stress relation, ex- 
pressed in these tilted axis, can be written in the following 
matrix form: 

U» = -ft & o U 2 , (i) 

W/ V is/Viz/ 
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Because the matrix involved in this relation must be sym- 
metric, the Young's moduli and Poisson ratios are not in- 
dependent, but satisfy the relation E%/Ei = ^21/^12- The 
elastic free energy is well defined if 



E u E 2 ,G>0, 

1 - ^12^21 > 0. 



(2) 
(3) 



The isotropic case is recovered for E\ = E2 = E, v\ 2 = 
v 21 = v and G = ^iqby- 

In reference Otto et al. have analytically com- 
puted the stress tensor components for such an anisotropic 
elastic material in the case of a semi-infinite medium (z > 
0) loaded with a unitary point force at the origin x — z = 
0. They have shown that, for a given anisotropy angle r, 
the stress response profiles only depend on two quantities 
that involve a combination of the elastic coefficients: 



t = E 2 /E 1 = V2l/V\2, 

1 „ / 1 ^12 ^21 
r = — E -2 \ — 

2 2 \G E x E 2 



(4) 
(5) 



It is notable that for other geometries, the stresses could 
depend on the E's, v's or G independently of t and r due 
to the effect of different boundary conditions, for example 
at the bottom of a finite thickness slab. By contrast, the 
corresponding familiar isotropic solution is completely in- 
dependent of E and v, although some (weak) variations of 
the stress profiles with v are found for an isotropic layer 
of finite thickness |22I48| . Depending on the values of t, 
and in particular on the sign of r, and of (r 2 — t), dif- 
ferent response shapes result, including single or double 
peaked responses, and symmetrical or skewed profiles (see 
numerous figures in 

Although it is natural to fix the angle r and to use the 
ratio E2/E1 as a control parameter for the amplitude of 
the anisotropy created by the shearing, it is rather difficult 
to have any kind of intuition about the i/'s and the shear 
modulus G, i.e. on r. Figure shows the locus of values 
of r and t which give a response profile for a semi-infinite 
medium with a deflection angle of (3 — 22°, the experimen- 
tal value. Interestingly, we see the tendency for a higher 
value of r to significantly broaden the response function. 
This would correspond, for example, to small values of the 
shear modulus G. Note also, that the width at half-height 
of the response is a quantity that is easily accessible ex- 
perimentally, and hence provides a good method to test 
this type of model. 



2.2 Elastic response in 3D 

It is not difficult to generalize this approach to 3D. How- 
ever, in this case, two axes (directions 2 and 3) orthogonal 
to direction 1 must be specified. We keep the notation x 
for the horizontal axis along which the shearing is applied, 
and the stress response measured. Direction 2 is in the ver- 
tical plane {z,x), and direction 3 (y) is perpendicular to 
these two directions. 
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Fig. 6. Values of r and t which give a response profile for a 
semi-infinite medium with a tilt angle of /3 = 22°. Inset: Three 
particular such stress profiles: r = 0.15 and t — 0.07 (A, bold 
line), r = 1.0 and t = 0.24 (B, thin line), and r = 2.5 and 
t = 0.62 (C, dashed line). The response shape is wider for 
larger r. 

The matrix involved in the 3D equivalent of the strain- 
stress relations Q has a similar structure to the 2D case, 
but with three Young moduli Ei, three shear moduli Gi 
and six Poisson ratios i/y (i and j must be different). 
Again, the symmetry of the matrix gives three relations of 
the form v^jEi = i/ji/Ej. In the final analysis, there are 
9 independent coefficients besides the anisotropy angle r. 
Stability now requires that 



Ei,Gi > 0, 
1 - fijVji > 0, 

^21^32^13 > 0. 



(6) 
(7) 

(8) 



Unfortunately, in 3D, we do not yet have closed form an- 
alytic expressions for the stress response profiles compa- 
rable to those given in PHj. To compute profiles in 3D, 
we make use of a finite element free-ware code developed 
in soil mechanics called CASTEM j3H|. We typically chose 
a grid with 40 x 40 x 12 cells in the x, y and z direction 
respectively. In this subsection, we show that the trends 
suggested by the experimental data are reproduced rather 
well in a model that allows for a modification of the con- 
stitutive elastic relations with respect to the symmetry of 
the external deformation applied to the system. 

Here, we analyze the previous 3D experiments on the 
shear box and on the avalanche preparation. We empha- 
size that we do not make a detailed fit of the experimental 
data to a general anisotropic elasticity model. Using the 
numerical code, we only explore the possibilities of or- 
thotropic elasticity on a 'Simple Elastic Model' that we 
call SEM in the sequel. Again, we do not seek a com- 
plete exploration of the parameter space and thus, we 
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Fig. 7. 3D equivalent of figurc|S| values of u and t which give a 
response profile with a tilt angle of /3 = 8° . Note however that 
these 3D calculations have been performed numerically with 
CASTEM on a layer of finite thickness. Inset: Three particular 
such stress profiles: u = 1.05 and t = 0.43 (A, bold line), 
u = 2.14 and t = 0.67 (B, thin line), and u = 3.96 and t = 0.8 
(C, dashed line). The response shape is wider for larger u. 



to = 0.67. The curve corresponding to (uq, to) is presented 
in figure 01 

The avalanche preparation - Now, we seek an inter- 
pretation of the response function for the granular slab 
prepared by successive avalanches. First we note that dur- 
ing the avalanching process, the flowing layer experiences 
a strong shear along the avalanche angle ip. An interest- 
ing outcome of the shear box experiment is that for all 
finite shear angle values <fi that we tested, the shape of 
the response function as well as tilt angle (3 did not vary 
significantly. We argue that each deposited layer retains 
a memory of the shear due to the avalanching process, 
which induces an anisotropic texture. In first approxima- 
tion, we may assume that, the avalanche acts like a shear 
box whose bottom is tilted at a angle ip ~ 30° with respect 
to the horizontal direction. We propose to use the param- 
eters (itn, tg) which compares reasonably to the shear box 
data to interpretate the avalanche deposition experiment. 
Of course now, we need to tilt the orthotropic angle r from 
45° to a value r = 45° — <p ~ 15°. In figure we present 
the results of CASTEM calculations on the SEM using these 
parameters and indeed we find an angle j3 = 3.0° ± 0.5°, 
close to the experimental value f3 = 3.5° ± 1°. 



3 Back to the sandpile pressure profile 



only vary three parameters that we believe are crucial 
namely, the stiffness of the direction of anisotropy E% , the 
shear modulus G and r the direction of the orthotropic 
axis with respect to the vertical direction. More precisely, 
we consider the same value of G for all directions, i.e. 
G\ = Gi = G3 = G. The average Young's modulus is 
taken constant: E = (E 1 + E 2 + E 3 )/3 = 150MPa. We 
also take E 2 = E3. Finally, the three Poisson coefficients 
^12, V13 and ^23 are also kept constant and bear the val- 
ues v = 0.3. We checked that the value of this parame- 
ter is not very sensitive. The remaining are such that 
the strain-stress matrix is symmetrical (see above) . In the 
spirit of the 2D analytical results, we use two dimension- 
less parameters to present our data, i.e. the stiffness ratio 
t = E2/E1, and the shear ratio u = E2/G. u is analogous 
to the parameter r in 2D - see equation (J5J. 

The shear cell - An essential result of the shear box 
experiment was that the tilt angle of the stress response 
profiles is around 8° (figure 0J. As suggested in the 2D 
space, it is natural that the orthotropic direction is tilted 
with an angle r = 45°. In figure we thus present the 
curve t(u) corresponding to a tilt angle (3 = 8° and for 
r = 45° as computed by CASTEM using the SEM. We see 
that similarly to the 2D situation, moving on this curve 
from lower to higher values of u significantly enlarges the 
width of the response function (see inset of figure 0). This 
is indeed an experimental observable that can be used 
to discriminate in parameter space. In the framework of 
the SEM, good agreement with the experimental data on 
the shear box can be obtained for values uq — 2.14 and 



Several experimental determinations of the stress distribu- 
tion below a sand pile prepared from a point source were 
obtained under various experimental conditions. First, we 
note that all stress data can be rescaled so that it is pos- 
sible to compare the prediction of the elastic model with 
experiments made with different pile sizes. For a conical 
pile of height h, radius at the base R = h/ tan ip, and pile 
slope angle tp, the relevant parameters are the rescaled 
vertical stress o~ zz /pgh and the rescaled horizontal posi- 
tion £ — r/R (p is the density of the granular packing and 
g the gravity constant). 

In figure |H1 we show the experimental data obtained 
by three groups, Smid and Novosad 0, Brockbank et al. 
[3] and Vanel et al. 0]. The underlying idea behind the 
use of an anisotropic elasticity is to include the frozen 
texture caused by shearing that ensues during avalanche 
deposition. From the response function experiments un- 
der shearing (shear box and avalanche) displayed before, 
it is natural to propose a modelling of the sand pile pre- 
pared by successive avalanche using an orthotropic elastic 
model like SEM. Of course, we see here the limitations of 
such a modelling, since in the SEM version of the CASTEM 
computation we only have three parameters to vary. Fur- 
thermore, as discussed previously, in the case where the 
sand pile repose angle is (p = 30°, the orthotropic direc- 
tion should be fixed at r = 15°. We show in figure |S] the 
results of a CASTEM computation for a conical sand pile 
using the 'optimal' parameters (urj,io) that yielded a rea- 
sonable agreement to the shear box and the avalanche de- 
position experiments. Note that we also checked that the 
stress state is always below the Coulomb criterion (with 
an internal friction angle of 30°), i.e. that failure never 
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— t=0.67, u=2.14 

— t=0.67, 
• h=20cm 
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Fig. 8. Normal stress below a conical pile. Filled symbols are 
the data of Smid and Novosad the stars are from Brockbank 
et al. [3], and the empty circles show Vanel et al. measurements 
The corresponding pile heights h are indicated in legend. 
The repose angles of the material used by these authors are 
respectively ip — 33°, cp = 31° and ip = 30°. For simplicity we 
have set ip — 30° in all our calculations. The solid line is the 
prediction of the SEM in the conical geometry with t = 0.67 
and u = 2.14 when r is set to 15°. The dashed line is the same 
except for u — > oo, which seems to fit better the data. 



occurs - except marginally at the very free surface as all 
stress components vanish. 

We observe the presence of a dip that bears essentially 
the correct qualitative features. Nevertheless, on a quan- 
titative level, the values of the depth of the dip as well as 
the amplitude of the stress maximum seem too large. Fur- 
thermore, the value of £ ma x ~ 0.15 underestimates clearly 
the experimental data. In fact the data seem to be better 
fitted by the SEM with a very small shear modulus, i.e. 
u —* oo. The understanding of this observation remains 
an open question. 



4 Conclusion 

In this paper we present an analysis of several response 
function experiments under the scope of anisotropic elas- 
ticity. As an application, we investigate to which extent 
the pressure dip observed at the bottom of a sand pile 
prepared by successive avalanches can be understood as a 
texturing effect induced by shearing along the avalanche 
directions. 

In the context of an elastic analysis, the orthotropic 
model used in this paper is relatively simple, since its main 
features are the presence of a stiff axis in the orthotropic 
direction and a shear modulus that can be varied indepen- 
dently. The orthotropy axis direction is a third parameter 
simply taken at an angle ir/4 with respect to the shear- 
ing direction. This choice corresponds to the direction of 



principal compression and is consistent with our recent ex- 
perimental findings in 2D . We also present new sets of 
experiments on a 3D sand packing in a shear box, where 
shear induced anisotropy is present. The second experi- 
ment is a sheared granular slab constructed by successive 
avalanches. The agreement with the orthotropic model is 
quantitatively correct, as the tilt of the response function 
can be reproduced by elastic modelling in both experi- 
ments using the same set of parameters. 

In a second series of analysis, we use the best fit pa- 
rameters obtained from comparison with response func- 
tion experiments to see whether such a relatively simple 
modelling is likely to explain the dip below a sand pile, as 
suggested initially by Savage . Indeed from a finite ele- 
ment calculation of a conical pile, we obtain a good quali- 
tative agreement with the available experimental data but 
we fail to obtain really quantitative agreement, as the dip 
amplitude and its width seem both underestimated by a 
factor of almost 30% in the central part of the pile. It is 
not yet clear whether the elasticity model we use is over- 
simplified since it contains only two parameters which are 
sensitive to shear, or whether the preparation procedure 
under avalanches is not accounted for correctly in the con- 
text of the 3D sand-pile. We note here the interesting sug- 
gestion of Jenkins 50 that the upward moving 'stoppage' 
waves produced when the avalanche hits the ground could 
modify the main compression axis so that its direction 
could be further away from the vertical. This effect would 
indeed enlarge the size of the dip. 

The tendency for deposition history or external action 
like shear or a biaxial compression to modify the constitu- 
tive structure of a material was noted in experiments by 
Oda et al. |51j and in numerical simulations by Radjai et 
al. |52j and is at present time still a very open and difficult 
question. This paper calls for more extensive systematic 
studies both experimentally and numerically (possibly an- 
alytically). The sand pile pressure dip, as interesting as 
it seems, appears a bit too complex to analyse for the 
moment. For a better understanding of this crucial issue, 
this work suggests a systematic use of response function 
techniques - stress responses as here, as well as displace- 
ment responses as in |53l54l55j - so as to extract several 
effective constitutive parameters of the material along a 
given stress-strain history. This is a promising systematic 
approach, and a possible alternative to sound propagation 
techniques, which can precisely identify internal structural 
changes due to external action a granular material |34j . 
This is the scope of ongoing projects in our laboratories. 
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